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We report first-principles studies of electronic transport in monolayers of Tour wires functionalized 
with different side groups. An analysis of the scattering states and transmission eigenchannels 
suggests that the functionalization does not strongly affect the resonances responsible for current 
flow through the monolayer. However, functionalization has a significant effect on the interactions 
within the monolayer, so that monolayers with NO 2 side groups exhibit local minima associated 
with twisted conformations of the molecules. We use our results to interpret observations of negative 
differential resistance and molecular memory in monolayers of NO2 functionalized molecules in terms 
of a twisting of the central ring induced by an applied bias potential. 
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An important goal in the study of molecular electron- 
ics is to identify molecules that can be combined to per- 
form logical functions 0. In particular, Reed and co- 
workers ||, ||, ||, [5| have studied the electrical proper- 
ties of a set of phenyl-ethylene oligomers (known as Tour 
wires (TW)||) functionalized with different side groups, 
and demonstrated that such molecules can show nega- 
tive differential resistance (NDR)||, || and a "molecu- 
lar memory" ^|, ^ effect, in which case molecules can be 
switched from a low conductance state to a high con- 
ductance state by application of a voltage pulse. A fun- 
damental understanding of the microscopic mechanisms 
governing the NDR and memory effects in these systems 
is still lacking. 

Complementary studies, by Weiss and co-workers [[fj, 
of isolated or small bundles of TWs embedded in a 
monolayer of spacer alkyl molecules, demonstrated spon- 
taneous conductance switching. While Reed and co- 
workers only observed switching effects for TWs with 
NO2 sidegroups, Weiss et. al. reported spontaneous 
switching of all the molecules in their study but had 
limited success in inducing switching by using a volt- 
age pulse [Q. Theoretical studies of properties of iso- 
lated TWs have emphasized the role of charging of the 
molecule and subsequent localization/delocalization of 
molecular orbitals ||, ^| or bias-induced alignment of 
molecular orbitals on the first and last phenyl rings of 
twisted TWs[|o|. 

In this paper we use the TranSiesta package jn| to in- 
vestigate the effect of different sidegroups on both the 
electrical and structural properties of TWs in self assem- 
bled monolayers, and in this way gain new insight into the 
microscopic origin of the observed NDR and memory ef- 
fects. TranSIESTA is a non-equilibrium Greens function 
based electronic transport program founded on density 
functional theory. It treats the entire system selfconsis- 
tently under finite bias conditions, and the only major ap- 
proximation in the method is the choice of exchange cor- 
relation(XC) functional. [O We have previously shown 
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FIG. 1: Geometry of monolayers A, B, and C connected 
with two Au (111) surfaces. Color codes: C (green), H (white), 
O(red), N(blue), S (yellow), Au(gold). 



that electrical properties of different metallic wires calcu- 
lated using the method are in excellent agreement with 
experimental data]l3[. 

The nanopore structures in which NDR and memory 
effects have been observed are formed by evaporating the 
molecules onto a small Au surface inside a SiaN4 struc- 
ture, and then evaporating another Au electrode onto 
the molecular monolayer [||. Since there is no detailed 
information about the structure of the monolayer, our 
theoretical analysis of the problem proceeds by forming 
an idealized model of the electrode/monolayer /electrode 
system and comparing electron transport in three dif- 
ferent molecules]!^]. We use an Au(lll) surface in a 
3x3 unit cell for the electrodes and assume the TWs 
are chemisorbed to the surfaces through strong thiol 
bonds [H, as illustrated in Fig. 1. monolayer A consists 
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FIG. 2: (a) Zero-bias tranmssion T(E, Vb = OV) versus inci- 
dent E for monolayers A, B, and C. HOMO and LUMO reso- 
nances are indicated, (b) I-Vb characteristics for monolayers 

A, B, and C (c) Isosurface of transmission eigenchannel at the 
HOMO resonance for monolayers A, B, and C (E — —0.84, 
—0.65, and —0.85 eV respectively) (d) Isosurface of transmis- 
sion eigenchannel at the LUMO resonance for monolayers A, 

B, and C (E = 1.56 , 1.59, and 1.14 eV respectively) 



of a bare TW with thiol end groups, while monolayers B 
and C have been functionalized with NH2 and NO2 side 
groups, respectively. 

We first investigate the effect of the sidegroups on 
the electrical properties of the monolayers. The current 
through the monolayers is determined by the quantum 
mechanical probability for electrons to tunnel from one 
electrode to the other, and is calculated using the Lan- 
dauer formula^: I = J*** T(E,V b )dE, where fi L/R = 
±eVb and T(E, Vb) is the transmission probability for 
electrons incident at an energy E through a device under 
a potential bias Vb (n.b Vb > corresponds to hole injec- 
tion from the right electrode). The general shape of the 
zero-bias tranmsission spectra, shown as T(E, Vb = OV) 

in Fig. 2a, is similar for all three systems. The zero bias 
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conductance G = ^T(E = /i// r , Vb = 0), given by 2.0 /iS, 
2.3 /iS, 1.9 /iS, for monolayers A, B, and C, is dominated 
by the tail of a broad highest-occupied-molecular-orbital 
(HOMO) resonance (at ^HOMO = -°- 84 > -0.65, and 
—0.85 eV for monolayers A, B, and C). The narrower 
lowest-unoccupied-molecular-orbital (LUMO) resonance 

(at ^LUMO = L56 ' L59 ' and L14 eV for A ' B ' and 
C respectively) contributes less to the conductance and 

thus holes are the dominant carriers. The NH2 group 
is electron donating, so that ^HOMO b ~ ^HOMO a ~ 
200 meV while the NO2 group is electron accepting so 

^LUMO,c - E LVMO,a ~ - 400 meV - 

To calculate the I-Vb spectrum we performed self- 
consistent calculations for biases in the range -2.8 V to 
2.8 V in steps of 0.2 V. We note that the charge on the 



molecule is not fixed ]ll[ and adjusts itself to minimize the 
free energy |l7j as the left /right electrochemical poten- 
tials are changed. For all three systems we find that the 
molecules are close to charge neutrality, and the charge 
on the molecule Qm , as determined by a Mullliken pop- 
ulation analysis, change by less than 0.05 e as the bias 
voltage Vb is varied. The I-Vb spectra, shown in Fig. 
2b, are very similar, increasing slowly at first but then 
increasing rapidly around 2V where the resonances come 
into alignment with the bias window. The main effect 
of Vb on T(E, Vb) is to sample more and more of the 
resonance while resonant peak height and position vary 
slowly with T4|@- 

When the molecules interact with the Au(lll) 
electrodes, the molecular levels broaden into a 
continuum [[HI. The eigenstates of the whole 
metal/monolayer /metal system consist of scattering 
states ]l9[ [20| , which are molecular orbital-like in the 
molecule, and Bloch wave-like in the metal slabs. If an 
orbital is delocalized across the molecule, an electron 
that enters the molecule at the energy of the orbital 
has a high probability of reaching the other end and 
thus there is a corresponding peak in the transmission 
probability T{E, V&), as illustrated in Fig. 2a. By 
calculating the continuum eigenstates at the resonance 
energies (^HOMO' ^LUMO)' tne or bitals that are 
responsible for current flow through the molecule can 
be analyzed. The HOMO resonances, illustrated in Fig. 
2c, resemble the HOMO of the isolated molecules while 
the narrower LUMO resonances, illustrated in Fig. 2d, 
resemble the LUMO of the isolated molecules. There 
are minor differences between the orbitals: the NH2 side 
group participates in the HOMO resonance but not in 
the LUMO resonance, while the NO2 group participates 
in the LUMO resonance and not the HOMO resonance, 
consistent with their respective donating/accepting 
characters. It is, however, clear from both the shape of 
the transmission curves in Fig. 2a and the orbitals in 
Figs. 2c and 2d that the resonant transmission peaks in 
monolayers A, B, and C are related to the delocalized 
nature of the 7r-orbitals of the bare TWs and is not 
strongly affected by the functionalization. 

The energetics of the monolayers are, however, 
strongly affected by the functionalization. The triple 
bond between the phenyl rings is rotationally symmetric 
and thus constitutes an easy axis about which the mid- 
dle ring can rotate. The differences in total energy as a 
function of the rotation angle of the middle ring, for 
the molecules arranged in the 3x3 unit cell of Au(lll) 
are shown in Fig. 3a. The intermolecular interaction 
energy is strongly dependent on the functionalization: 
there is a cost of 60 meV (180 meV) to rotate the ring 
90°in monolayer A (B), and the flat molecule is clearly 
favoured. Monolayer C shows a quite different behaviour, 
it acquires a local minimum at = 60° and 120°. To un- 
derstand the origin of these minima we have plotted the 
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FIG. 3: (a) Energy versus rotation angle for Molecules A, 
B, and C in Au(lll) 3x3 unit cell. The energy is calculated 
within the Perdew, Burke and Ernzerhof (PBE) approxima- 
tion for the exchange- correlation functional[^. (b) Contour- 
plot of the effective potential between TW's with NO 2 side- 
groups. Note the bond formation between the O atom and 
the H atom on the neighboring TW. 

effective potential within the monolayer in Fig. 3b. The 
effective potential reveals the formation of a hydrogen 
bond between the NO2 group and a hydrogen atom on 
the neighboring TW. Similar hydrogen bond formation 
involving NO2 groups has also been observed in other 
molecular layers [^2j. From Fig. 3b it is also clear that 
the minima at 60°and 120°can be attributed to bond 
formation with a new neighboring molecule, and these 
particular angles arise from the symmetry of the mono- 
layer. 

There is a barrier for rotation of the middle ring from 
= to = 60°of ~ 150 meV, while the barrier from 
60°to 0°is ~ 65 meV We postulate that, as a voltage 
pulse is applied, these barriers may change, and a tran- 
sition to a new molecular conformation with the middle 
ring twisted may be induced. The local minima only ex- 
ist when the molecules are functionalized with an NO2 
side group and thus such transitions only exist when the 



molecules are functionalized with an NO2 side group [[|. 

Again because of the lack of detailed information about 
the monolayer, we have selected a candidate conforma- 
tion, which we designate monolayer C@60, formed by 
twisting the middle ring of molecules in monolayer C 60°, 
for possible stabilization under applied bias. The trans- 
mission spectra for monolayer C@60 is shown in Fig. 4a 
and its I- V& characteristics in Fig. 4d. The conductance 
G ~ .12/iS is roughly 16 times smaller than monolayer C 
while the current at V& = 2V is « 5 times smaller. To 
lowest order, G is proportional to the product of the ma- 
trix elements between 7r-orbitals on neighbouring phenyl 
rings so that a rotation of the middle ring will reduce G 
by a factor ~ cos 4 (#) = |, in rough agreement with the 
calculated values. 

The transmission eigenstates for electrons travel- 
ling from left to right at ^HOMO = — 0-85 eV and 
^LUMO = °' 99 eV are illustrated in Fig. 4b. The 
HOMO resonance is localized on the phenyl ring near- 
est the left electrode while the LUMO state corresponds 
to an orbital of the middle ring, weakly coupled to the 
first and last rings. 

To investigate whether monolayer C@60 could be sta- 
bilized with respect to monolayer C as the bias voltage is 
increased, we show in Fig. 4b the Kohn-Sham total en- 
ergy of monolayers C and C@60. At V*> = 0, monolayer C 
is lower in energy than monolayer C@60. However, when 
the bias is increased above 2 V, monolayer C@60 becomes 
lower in energy, leading to an expected conformational 
change. While there are issues surrounding the evalua- 
tion of total energies for such systems ^3[ , Fig. 4b 
serves to illustrate that different conformations may be 
stablized with respect to each other by an external bias. 
In this case, the calculation suggests a transformation 
of monolayer C to monolayer C@60 at around 2V. The 
exact value of the switching voltage will depend on the 
coverage, detailed structure and size of the monolayer, 
as well as how the potential drops across the molecules, 
which can be affected by electrode coupling ]l8[. 

The local minima in Fig. 3a are stabilized by the NO2 
sidegroup, which correlates well with that the molecular 
memory effect is only observed in monolayers with NO2 
side groups [Q. Furthermore, the energy barrier for ro- 
tation of 65 — 150 meV derived from our crude model 
is in rough agreement with the 80 meV barrier height 
extracted from bit retention times [Q, thus conforma- 
tional changes stabilized by intermolecular interactions 
is a good candidate for explaining such phenomena. 

Because the minima emerge from interactions within 
the monolayers, the above effects may only be observed in 
nanopore experiments ||, ||, |[ and not in experiments 
of single molecules. In the STM experiments of Weiss and 
co-workers on single or small bundles of molecules, the 
conductance switching may be related to different confor- 
mations formed by steric interactions between the TWs 
and the insulating alkyl molecules. This could explain 
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FIG. 4: (a) Zero-bias transmission T(E, Vb — 0) for mono- 
layer C@60. (b) Transmission eigenchannels corresponding to 
HOMO and LUMO resonances of C@60. (c) Total energy of 
monolayer C(blue) and monolayer C@60 (magenta) as a func- 
tion of bias potential, (d) /- Vb characteristics of monolayer 
C(blue) and monolayer C@60 (magenta). 



why the observed switching behaviour was independent 
of the functionalization of the molecules. 

In conclusion, we find that functionalization of TWs 
has a stronger effect on the energetics of the monolay- 
ers than on the orbitals responsible for current transport 
and a better understanding of the intermolecular interac- 
tions in such monolayers could hopefully be exploited in 
order to design molecular electronic devices with specific 
properties. 
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